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ON THE WEAK-COUPLING LIMIT FOR BOSONS AND FERMIONS 
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Abstract. In this paper we consider a large system of Bosons or Fermions. We start with an 
initial datum which is compatible with the Bose-Einstein, respectively Fermi-Dirac, statistics. 
We let the system of interacting particles evolve in a weak-coupling regime. We show that, in 
the limit, and up to the second order in the potential, the perturbative expansion expressing 
the value of the one-particle Wigner function at time t, agrees with the analogous expansion 
for the solution to the Uehling-Uhlenbeck equation. 

This paper follows in spirit the companion work [2], where the authors investigated the 
weak-coupling limit for particles obeying the Maxwell-Boltzmann statistics; here, they proved 
a (much stronger) convergence result towards the solution of the Boltzmann equation. 


1. Introduction 

In 1933 Uehling and Uhlenbeck in Ref. [17] proposed the following kinetic eqnation, 
called U-U in the seqnel, for the time evolntion of the one-particle Wigner fnnction /(x, v, t) 
associated with a large system of weakly interacting Bosons or Fermions (see Ref. [18] for 
the dehnition of the Wigner fnnction). The U-U equation is 


dtf{x,v;t)+v-'\/j;f{x,v;t)= / dn* / dn' / dn'lF(n, n*|n', n') 


f'fUl + S 7 r ^ ef){l + - fMl + 87r^d/')(l + STT^d/:)), 


( 1 . 1 ) 


where we use the standard short-hand notation 


f = f{x,v,t), /* = /(x,n*;t), f'= f{x,v'-,t), /'=/(a;, n'; t). 

Here, W denotes the transition kernel 

lF(n,n*|n',n') = i [^(t' - v) +e^{v' -v^)f 

1 ( 1 - 2 ) 
S{v + v^ -v' -v'^) + vl -v'"^ 
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Finally, 

(^{k) = J dx (1.3) 

is the Fourier transform of the two-body interaction potential (/>, and 9 = ±1 for Bosons 
and Fermions respectively. 

Note that the factors do not appear in the original U-U equation in Ref. [17], 
because there, the distribution function is normalized in such a way that its integral on 
the velocity variable equals the space density times Stt^. At variance, in (1.1), / is just 
the standard Wigner function, whose integral on the velocities equals the space density. 
Let us mention that equation (1.1) actually is cubic (and not quartic) in the unknown /: 
apparent quartic terms have vanishing contribution, as shown by direct inspection. 

Eq. (1.1) constitutes a natural modihcation of the usual quantum Boltzmann equation, 
in order to take into account statistics. In particular, there is a i7-functional 


nf) = 


dv{f log f -9(1 + 871^9 f) log(l + 87r^9f )} 


driving the system to the Bose-Einstein and Fermi-Dirac equilibrium distribution: 


M{v) 


^(3fi+Pv'^/2 _ 


outside the Bose condensation region. Here f3 and (i denote the inverse temperature and 
the chemical potential respectively. 

Eq. (1.1) is largely studied (see for istance [1] and [12] for physical consideration, and 
[15], [7], [13], [14] ... for a more mathematically oriented analysis concerning the existence 
of solutions and asymptotic behavior), so that it is certainly of great relevance to derive 
this equation from the hrst principles, namely from the Schrodinger equation. 

As clearly explained by H. Spohn in [16], Eq (1.1) is indeed expected to hold in the so 
called weak-coupling limit, which consists in scaling space, time and the potential according 
to 

X —> ex, t ^ et, (/) —> (1-6) 

where £ is a small positive parameter. 

A slightly different limit, usually called van Hove limit, scales t and cf) as in (1.6) but 
leaves the microscopic space scale unchanged. Eq. (1.1) cannot be derived in the van Hove 
limit in general but, in case of translationally invariant states, we expect to achieve the 
homogeneous version of the U-U equation (for a large system). In fact Hugenholtz [11] 
proved formally that this happens. Later on Ho and Landau [10] proved that the homo¬ 
geneous U-U equation holds rigorously up to the second order expansion in the potential. 
These approaches, as well as the recent contribution by Erdos et al. [8] (where the quantum 
analog of the Boltzmann’s Stosszahlansatz is formulated), are based on the commutator 
expansion of the time evolution of the observables of the CCR and CAR algebras. 

In the present paper we approach the problem from a different viewpoint. We start 
from the time evolution of a particle quantum system in terms of the Wigner formalism. 
Here the statistics enters only through the choice of the admissible states we take as initial 
conditions. Such states, called quasi-free, must describe free Bosons and Fermions, so that 
they cannot have any other correlations but those arising from the statistics. Therefore 
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the first step is to characterize quasi-free states (see for example Ref. [4]) in terms of 
the Wigner functions. Then we apply the dynamics (in terms of the usual hierarchy) and 
represent the solution as a perturbative expansion. The truncation of this expansion up 
to the second order in the potential is shown to converge to the expansion associated to 
the U-U equation, up to the hrst order in the scattering cross section. 

In other words we recover the result in Ref. [10] with the following main differences. 
First we exploit the weak-coupling limit, so that we can deal with states which are not 
necessarily translationally invariant. Second, we work directly in terms of the Wigner 
formalism, in the same spirit of the Balescu book (see Ref. [1]). In doing so, we also 
follow a previous work [2] by the authors for the Maxwell-Boltzmann statistics. Hence the 
present work shows how the statistics can be handled in this formalism. Note in passing 
that the case of the Maxwell-Boltzmann statistics allows for a much stronger (but still 
partial) convergence result than the one presented here, see [2]. Note hnally that the 
present formalism also allows to handle the low-density limit, see [3], see also the last 
section of this text. 

It is also important to mention that a full rigorous derivation of the U-U equation (but 
also of the usual Boltzmann equation arising for the Maxwell-Boltzmann statistics) is still 
far beyond the present techniques and those of the previous references. 

The plan of the paper is the following. In the next section we describe the particle 
system. In Section 3 we establish the result. The rest of the paper is devoted to the 
proofs. 


2. The particle system. 

We consider a Quantum particle system in Let 

i3 = 0T2(M3)-:=0i3n, (2.1) 

n>0 n>0 

be the Fock space. A state of the system is a self-adjoint, positive trace class operator 
acting on Sj: 

a = an- (2-2) 

n>0 

We assume 

Tra = l. (2.3) 

The operator N, number of particles, is the multiplication by n on Sjn and hence 

(TV) = J^nTra^, (2.4) 

n>0 


where the left hand side is the average number of particles in the state a. If an{Xn', W) is 
the kernel of the Reduced Density Matrices (RDM) are dehned by: 


Pn(A^n) Ui) ^ ^ 


(n -|- m)! 


i>0 


mi 


^m) dZfyi. 


(2.5) 
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Here Xn = (ti, ..., Xn), Xi G denotes the n-particle configuration. Note that 
Tl’Pn / dZ^i 

( 2 - 6 ) 

= 22 — 1).. . (m — n + l)TTam = {N{N — 1).. .{N — n + 1)), 

m^n 

and hence the RDM are equivalent to the classical correlation functions. 

The Hamiltonian of the system is the self-adjoint operator acting on S) given by 

CXD 

H = (2.7) 

n=l 


where 

1 ” 

= “2 ^ (2.8) 
j — 1 

and the potential (p is a, smooth two-body interaction. Here, h as well as the mass of the 
particles are normalized to unity. 

Under these circumstances, the time evolved state is given by the usual 

a(t) =e-*^Ve*^*. (2.9) 

Now, quantum statistics is taken into account by suitable properties of the physically 
relevant states. Namely, for the Maxwell-Boltzmann (M-B) statistics we require symmetry 
of pn{xi,... ,Xn;yi ■■■ ,yn) in the exchange of particle names. For the Bose-Einstein (B-E) 
and Fermi-Dirac (F-D) statistics we require additionally 

Pn(Tl, . . . , Tn; Pl . . . , I/n) = . . . , X^, y^{l) ■ ■ ■ , 2/^(n)), (2-10) 

where tt E Vn is a permutation of n elements, and s(7r) = 0 if the permutation is even, 
s(7r) = 1 if it is odd. 

Alternatively, the quantum statistics is automatically taken into account by considering 
states on the algebra generated by the annihilation and creation operators a{x) and a^(a:) 
(with the commutation and anti-commutation relations according to the B-E and F-D 
statistics respectively). Then the RDM are defined as 

Tr [aa\xn).. .a\xi)a{yi).. .a{yn)] = Pn(Ti,...,Pi..., y^)- (2.11) 

However we do not use here the second quantization formalism. 

Given a state cr, we define the Wigner transform [18] by 

lUn(X^;K):=^^^y' + (2.12) 

Therefore the analogous of the classical correlation functions are the j-particle Wigner 
functions defined through 

F,iXj; V,) = Wj+^{X^,X^- U„ K). 

n>0 J J 


(2.13) 
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Note that the Fj’s are the Wigner transforms of the RDM pj, as one can easily check. 

Dne to the dynamics imposed by (2.9), it is a standard compntation to check that the 
Wigner f un ction evolves according to the Wigner-Lionville eqnation 

n 

dtWr, + = T^W^. (2.14) 

i=l 

As a conseqnence, the j-particle Wigner fnnctions Fj’s satisfy the associated hierarchy 


j 

dfFj + '^^Vi ■ VxiFj = TjFj + Cj+iFj^i, (2-15) 

i=l 

where Tj and Cj^i will be dehned below after Eq.(2.22), and the index j takes any valne 
between 1 and N. Eqnations (2.15) are analogons to the nsnal BBGKY hierarchy for the 
classical systems and are derived in a similar manner. Note that by the dehnition of the 
RDM the coefficient in front of Cj+i is one instead of — j. 


We now want to analyze (2.15) in the weak-conpling regime (1.6). Therefore, we set 




(2.16) 


where £ > 0 is a small parameter, and we scale the potential as well, by setting 

(j) (2.17) 

The resnlting, scaled, eqnation is 




i=l 


where 


(TJf’)(Xf,Vi) = 

0<k<£<j 


(2.18) 


(2.19) 


and the T^^’s are dehned as follows: if j = 1, we simply have Tf = 0; otherwise, 

/ h h \ 

tj (^Ti, ...,Xj]Vi,...,Vk- cr-,.. .,n^ + a-, 

On the other hand, the in (2.18) is compnted as: 


( 2 . 20 ) 




k=l 


( 2 . 21 ) 
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where 




fj+i (^Ti,..., Xj+i'Vi, ...,Vk-cr^,..., Vj+i + a|). 

Note that Tj and Cj^i are T? and C'j_|_i for e = 1. Last, we hx an initial condition seqnence 

{/o^,,}r=i (2.23) 

according to the qnantnm statistics, and perform the limit £ —0 in the resnlting system. 


h' 


( 2 . 22 ) 


Remark : Since 

J fo,i{^T'>j)dxdv — e^{N), (2.24) 

reqniring ||/o iI|li = 0{1), implies (N) = 0{e~^). In other words we are working in the 
Grand-canonical formalism and the density is antomatically rescaled. □ 

In the following we shall hx /q to be a given (independent of e) probability density /q. 
This means that its inverse Wigner transform 

P%x,y) = J dve^^'^ (2.25) 

namely the one-particle rescaled RDM, is a snperposition of WKB states. 

We now make assnmptions on the initial state to take into acconnt the statistics. For 
the M-B statistics a snitable initial seqnence can be chosen completely factorized, e.g. 

/o', = /o®/ (2-26) 

Snch a notion of statistical independence, which corresponds to a complete factorization 
of the RDM’s, is not compatible (bnt for the condensed Bose state) with the B-E and F- 
D statistics which exhibit intrinsic correlations even for non interacting particle systems. 
States describing free Bosons or Fermions are nsnally called qnasi-free and are dehned in 
terms of the RDM’s by the following formnla: 

j 

Pj{x\i ■ ■ ■ 1Xj^ yi,..., t/j^ ^ ^ 9 ^ ^ p(^Xi, yn{i)^■> (2.27) 

i=l 

for some positive dehnite operator p on L 2 (K^) with kernel p{x, y). We show in Appendix 
how to constrnct explicitly qnasi-free states for Bosons. 

From now on we assnme that the initial seqnence (2.23) for the rescaled problem (2.18) 
is given by the Wigner transform of a qnasi-free state (2.27) generated by p{x^ y) = p®(t, y) 
given by (2.25). As a conseqnence the initial seqnence for the hierarchy (2.18) is 

of the form 

f°(X„Vi)= Y. V,), 

neVj 


(2.28) 
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with 




1 


n 

k=l 


^iyk-Vk+iwk- 


xu—x^ 


dyv ■ ■ 

7ol 


{2'kY^ 

^ pfc+T^(fc) £ 


dyj / dwi - ■ ■ / dwj 


A {yk yn(k)): '^k 


(2.29) 


Eq. (2.29) follows from (2.27) and (2.25). 

We underline once more that, in the present approach, the dynamics is given by the 
hierarchy of equations (2.18) which are completely equivalent to the Schrodinger equation, 
while the statistics enters only in the structure of the initial state. 

In the weak-coupling limit £ —0, we expect that fj{t) converges to a factorized state 
(because the effects of statistics disappear in the macroscopic limit). On the more each 
factor should be solution to the U-U equation (the collisions being affected by the statistics 
because they involve microscopic scales). 


3. The main result. 

Let / = f{x,v,t) be a solution to the U-U equation and set /j( ■, ■ ,t) = 
Then the sequence {fj}^i satishes the following hierachy of equations: 

j 

(dt + ^ ~ Qj,j+^fj + ^ “k Qj,j+‘^fj+‘^- 

i=l 




(3.1) 


Here the Qj j^i contribution, a ’’two particles term” in the terminology used below, is 


= / dv'k / dvj+i / dVj+^W{vk,Vj+i\vk,Vj+^) 


k=i ' 


(3.2) 


[fi+i ti'+i) - fj+i(Xj,xi,;vi,..., 


and the Qj,j +2 contribution, a ’’three particles term”, is 


iQj , j + 2fj + 2 ) iXj , Vj ) = 871 ^ 9 ^ j dv'k j dvj+i j dv '^^^ W { vk , Vj + i \ v ' k , v ' j ^^) 

k=l 

/j+2(Wj, Xfc, Tfc; m,..., 4,... n'+i, Tfc) + 7+2 (Wj, Xfc, Xfc; m,..., 4,... n'+i, Tj+i) 


fj-\-2{Xj^ Xk Xk] t’l, . . . , Vj-\-\^ Vj^ fj-\-2{Xj Xk’i Xk] t’l, . . . , /. 


Also, (Xn, y) denotes the (n -|- l)-sequence (xi,..., x^ y). 


(3.3) 
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A formal solution to the hierarchy (3.1) is given by the following series expansion: 


m-Y, Y 


l>0 “1- -“n 


ft ptl 

dti / 

Jo 


f'tn-l 


dtr 


'0 


— oci=l,2 

S(t S{tl - t2)Qj+ai >i+cn+Q2‘^(^2 ^ 3 ) • • • 

O) S'O 1 -tan) 

Vj+aiH-|-an-l,i+aiH- hctn'^K^n) J Q : 

where S{t) denotes the fream stream operator, namely, 


(3.4) 


(3.5) 


As for the solution to the ^-dependent hierarchy (2.18), we can also expand /| in the 
similar way, at least at the formal level. This gives 


fm=Y Y 


n>0 T'l...T'n 
— 74 = 0,1 


ft fti 

dti / dt2- 
Jo 


— 1 


dtr 


'0 


S{t - -ts)... 


(3.6) 


P! 


j+7iH-l-7n-i,j+7iH-l-7i 




j+7lH-|-7n ’ 


where /• is an initial quasi-free state given by (2.29), and 


^7,.+i = ‘=“’0+.. 


(3.7) 


We are not able to show the convergence of fj{t) to fj{t) in the limit £ —0 even for 
short times. However we are going to show that the two series agree up to the second 
order in the potential. Namely, dehning the second order contributions 


g{t):=Sit)fo+ [ dtiSit-ti)Qi, 2 Siti)f^^+ [ dtiSit-ti)Qi,3Siti)f^^ (3.8) 
Jo Jo 

associated with /j(t), and 


ft 


g^{t) := S{t)fo + e-- dti S{t - ti)CIS{ti)f^ 

Jo 

ft pt 

,-4 


0 


+ £' 


+ £ 


-7 


dh / dt2 S{t-ti)CIS{h-t2)T^S{t2)fl 
Jo 

rti 

dti / dt2S(t-ti)ciS(ti-t2)cis(t2)fl 


(3.9) 


associated with /|(t), we rigorously prove below the convergence of g^{t) to g{t), under 
suitable assumptions on the data of the problem. 

Assumptions: We require (f) to be real and even, so that (ft is real. In particular 


m = ct>{-k) = ct>{-k). 
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This the most important assumption we need in the analysis. Besides, we shall need to 
deal with ’’smooth” data. Quantitatively, we assume the following regularity: 

(i + kl)" 

I/3|<2 

for a sufficiently large a, and 
fo{x,v) e Li, 

(i + iei + i>?ir E E i(of + o,';)/oK.>?)ieii. 

0<|/3|<2 0<|7|<2 

for a sufficiently large a as well. In (3.10), (3 and 7 denote multi-indices, and 
denote derivatives with respect to the variables ^ and rj. Note that throughout this paper 
we use the following normalization for the Fourier transform: 


JR" 

m = f dh /(A)e‘'-. 

1^'^; JR" 


(3.11) 


□ 


Our main result is the 

Theorem. Under the above assumptions, we have 

lim rj, t) = p, t), for any t > 0 and any (^, p) e M®. 
£^0 


Remark : In the above statement (and the proofs given below), we found convenient to treat 
the terms in (3.8) and (3.9) in terms of their Fourier transforms, for which the convergence 
arises more naturally. However, we would like to stress that in the companion paper [2], a 
stronger, but analogous, result is formulated in terms of the pointwise convergence in the 
X — V space, hence without going to the Fourier space. □ 

Before entering the details of the proof we hrst analyze all the contributions in the right 
hand side of (3.9). 

The two-particle terms are (we skip the unessential operator dti S{t — ti)) 

(3.12) 

where the permutation tt may take the two values tt = ( 1 , 2 ) or tt = ( 2 , 1 ), together with 

r dt2 CISih - t2)TfSit2)ff, (3.13) 

Jo 

with TT taking the values tt = (1, 2) or tt = (2,1). There are four such terms. 
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The three-particle terms are twelve, namely: 

Wls ■■= r *2 - i2)CpS(i2)/J, (3,14) 

Jo 

and ^ 

VI, := f ' dt2 cps(ti - t2)Ci,S(h)f^, (3.15) 

Jo 

with Ti E Vs, the set of the permntations of three objects, whose cardinality is six. 

Note that all the above terms are fnntions of {x,v) (and ti of conrse). 

For fnrther convenience, and in view of the proof of onr main result, we readily express 
all these contributions in terms of their Fourier transforms. 

We start with the following obvious three formulae for the basic operators S{t), T 2 , and 
C 2 (see (3.5), (2.19)-(2.20), and (2.21)-(2.22), respectively): 


T2fi^l,^2-,Vl:V2) = -i (J f dh 

a=±l 


m 

(27r)^ 




f {Cl-^,C2 + ^■,Vl,V2 ] , 


Clf{^]V)=T2f{^,0]V,0) = -i ^ / dh 

J ( 2 ^) 


Hh) 


h h 


e —0 , 


£ £ 


These relations give in (3.12) through (3.15): 


= -i 


_ 7 
8 2 


,-4 


^ cr I dh^{h) e *" ^ + fi(^ - j), j 


r = ±l 


£ £ 


(3.16) 


^2^,si^,v) ^ I dt2 I dhi I dh2${hi)^{h2) 


(Ti,(T 2=±1 


Tn [ c ^ (h ^ h \ ^ (U ^ h \ _L+t ^ 2/12 tihi t2h2 

/2 U - -(/ii + / 12 ), - [hi +h2)-,v + -,- 

\ £ £ £ £ 

^ Y] [ dt2 / dhi [ dh2 J(/li) J(/l2) 

Jo J J 

t-iT ff l/i, , , N ^2 , , . tihi+t2h2 tihi t2h2 

fs U - -Wi +/12), —, —+ -,-,-, ,, 

\ £ £ £ see 

V3,e(^i^) = - TA-Te Y ^1^2 [ dt2 [ dhi [ dh2 ^[hi) ^[h2) 

Jo J J 

- / hi hi-h2 h2 , ^ . tihi tihi-t2h2 ^ 2^12 

■^3 + -—,-;- 

V £ £ £ £ £ £ 


(3.17) 


(3.18) 


(3.19) 



ON THE WEAK-COUPLING LIMIT FOR BOSONS AND FERMIONS 11 

Starting form those expressions, the plan of the proof is the following. In Section 4 
we evalnate the two particle terms iSJg and We prove that they converge towards 

the associated two particles terms in the U-U eqnation. In Section 5 we deal with the 
three-particle terms associated to the permntations tt with a hxed element. Those are 
shown to converge towards the associated three particles terms in the U-U eqnation, while 
contribnting by the qnantity 4>{v' — vY + the transition kernel W (see (1.2)). 

Finally in Section 6 we treat the three particle terms relative to cyclic permntations. We 
recover in this way the missing contribntion to the transition kernel, namely the cross term 
64>{y' — v) 4>{v' — n*). 

For sake of simplicity we shall carry ont the compntations for Bosons (^ = 1), being 
clear that the Fermionic case is jnst the same with snitable changes of sign. 


4. Two-particle terms. 


We introdnce the partial Fonrier transform 


- J dvi- ■ ■ J dvj e fJ(xj,...,Xj-,Vi,...,Vj). 

(4.1) 


As a conseqnence of (2.29) we have 


II/o 


k=l 


Xk T X^Q^'^ lyfc Vn(k) ^k ^TT(k) Vk T VTT(k) 


— e 


+ 


2 4 ’ £ 2 

In particnlar, we have the obvions 

ff ’^\ xi , X 2 ; m , m ) = 7o(ti,?7i)7o(t2,i72), 

together with 

7(2,1). X jfxi+X2 e X 1 -X 2 vi+r]2 

H (ti,T 2 ; 771 ,^ 2 ) =/o I - ^ - |(hi -V 2 ); -^^— 

r f xi+ X2 s. X xi - X2 + m 

h (^— + jhi - 02 ); +^— + 

Hence, npon now performing the complete Fonrier transform, we obtain, 

together with 

ff-"\^ub-,riuV2) =e^ f dyi 


fo\yi- |(<7i - m)\ 
fo{yi + |(i7i -172); +!/2 + AilM. 


(4.2) 


(4.3) 


(4.4) 


(4.6) 


(4.6) 


2 
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We are now in position to analyse the term iSJ^ for tt = (1,2) and tt = (2,1). First, 
nsing the identity 

cr = ibl 

we get the the explicit expression: 


2 ’ 


^ yd/i J(/i) sin f 2 - 7 ’ 7 ;^+ (^ - 7 )^ 1 ’ 7 ^ 1 ) • (4-7) 

In the case of a change of variable h ^ eh then gives, nsing (4.5), the valne 


= “-^77^3 J dh(j){eh) sin ^ 

Therefore, we may estimate 


2e 2 


h ■ rj 


foi^ - h;rj + - h)ti)fo{h, hti). (4.8) 


i#T’i < 


< Cy/e\\^\\L^ I dh\h\ (|?7+ (^ - h)ti\ + |^ - h|ti) \fo{h;hti)\ |/o|(^ - h; ry + (^ - h)ti] 

(. 


snp|?7||/o|(^;?7) / t^^snp |^||/o|(^;?7)+ ^isup |^||/o|(^;?7) / snp |/o(^;?7)| 

\^,v J V J V 


< C^/esnp ({\C\ + |?7|) |/o|(^;h)j j d^ snp (^(|^| + |?7|) |/o|(^;h)j , 

and the corresponding contribntion vanishes with £. 

In the case of on the other hand, eqnations (4.7) and (4.6) give 


(4.9) 


4TK.») = 




d ' ^-i(y-^ + Ly2)-Fih-y2 


f e 

2h' 

\y^~ ~A 

h + - — 

\ 4 

eti 

f e 

2h' 

("‘ + 4 



V + Ch 
2 

V + ^ti 


;-2/2 + 
;2/2 + 


- j j j I e 

7 f ,d , V + ^ti 

fo\yi+ 2^1; -y^ ^ 2 — 


-i{yi+e^)-i^ih-y2 


yi - 7^1 W2 H-^— 


+ o(yi). 


(4.10) 

By the parity of (/>, the hrst term in the right hand side is vanishing: Indeed, it is anti¬ 
symmetric in the exchange h ^ —h and 1/2 —y 2 - Note that the mechanism that makes 

the dominant, 0(£“^/^), contribntion of S^d\ vanish in the limit, is very different from 

the one involved in the vanishing of iS;) J : here, antisymmetry plays a crncial role. This 
aspect will play an even more important, and more intricate, role in the next two sections. 
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There remains to prove that the 0(^/e) term in (4.10) indeed has the claimed size. It can 
be written as 


~ $(h) fo -y2^ fo +y2^ 

^ih-y2+i^-h-(2^l—^) (“2/2•?+(?—2?i^ | ^ ^ 

It may be estimated by 


Ce^ / dh\(f)\{h) / di/2C^6l/ol Ui; 


V + ^h 


+ 1 / 2 ] l/ol ( ^ — ^ 1 ; 


?7 + ^ti 


-2/2 


16 


?7 + 6i 


+ 2/2 


+ ie-6l 


?7 + 6 i 


- 2/2 


Therefore the term vanishes as well. 

As a conclusion, all terms S 2 g, which are the ones that are linear in (/>, vanish in the 
limit £ —0. 


We now pass to the evaluation of the terms T^g- 

The contribution has been already considered in Ref. [2]. However, for sake of 

completeness, we analyze this term in the present context as well. Using (4.5) in (3.17), 
and performing the change of variables: 

= = (4.11) 

S e 


we arrive at 
db2) 


r/b’K.’l) =-jAe ^ 1^2 ! ds I dh I dk${h)${-h +ek) 
(^’">+,.+±1 Jo J J 

Q-ilvh{ )+a- 2 h'^s] ^-ie^k-r] g-iso- 2 -I• [^s-2sfc] 

foi^ - k-,r] + tiC + hs - kti) fo{k; kti - hs). 


(4.12) 


This term converges formally to 

--{1,2).. . 1 


^ 2 ' ’^(6^7) = 


(2t) 


E 


'* + 00 


aia2 


ds dh dk \(f>{h)\^ 


(Ti,(T2=±1 


(4.13) 




To justify the limit we split the integration in ds over the two intervals [0,1], [1, +cx)]. In 
the hrst interval we bound the integrand by 


II<^IIl„ II/oIIi„ l<^(A)|sup l/o(£t?)l. 


(4.14) 
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which is a Li{dkdh) function for any s G [0,1]. In the second part of the integration 
domain, after the change of variables h {kti — hs), we bound the integrand by 

llJllLll/olU»l/o(*.'l)l4. (4.15) 

o 

which is a Li{dk dh ds) function on X X (1, +CX)). The claimed convergence in (4.13) 
is then consequence of the Dominated Convergence Theorem. It holds uniformly in rj. 

We now evaluate Inserting (4.6) in (3.17), and rescaling time ti — t 2 = ss, the 

resulting expression is: 


Y] ^1^2 / ds dhi / dh2 / dyi / dy2 ^(hi) ^{h2) 
a,,a2=±l Jo J J J J 

-i^h2-{ri+es^-2shi) + f ^ 2 ) „i{hi+h2)-y2 


fo\yi-^ 


i'l+i 




V + Ch 


hiti + h 2 {ti - es) 


hiti + h 2 {ti - es) 


; -2/2 + 
;2/2 + 


2 

V + ^ti 


Now the formal limit is: 


(4.16) 




1 


(27r; 


6 


E 


r*H-oo 


CTia2 


ds / dhi / dh 2 / dyi / dy 2 (t>{hi) (f){h 2 ) 


'0 


ai ,cr2 = ±l 

^-i^h2-{v-2shi) p,i{hi+h2)-y2 


fo ( 2/1 + ^1 1 ; -2/2 + 


e ^ e 

V + ^ti 


hi + /i2 , ^ + ^^1 

2/1 — ti --; 2/2 H-— 


(4.17) 

To justify the limit we have to show the uniform (with respect to e) integrability of the 
integrand in the right hand side of (4.16). To outline the decay with respect to the s 
variable we observe 


^icF2h\-h2S _ 


1 


s‘^\h2 




(4.18) 


and then proceed with the natural integration by parts with respect to hi in (4.16) (Recall 
that l/|h 2 p is integrable close to the origin in dimension d = 3). Splitting the integral in ds 
as before, we may apply the Dominated Convergence Theorem, upon using the smoothness 
of (f) and /o, thus justifying the above formal limit. 


Our last task is to interpret the result we have obtained, in terms of the U-U equation. 
To do so, we hrst go back to the original variables, expressing and as functions 

of (a:,u). A straightforward computation yields, on the one hand. 


^2 


( 1 , 2 ) 


{x,v) 


( 27 r) 


’Y) cricr2 J J *^^1 J *^^2 S ^ — Vi — h 


(72 — Ui 


(Tl ,cr2=±l 


A+ {-h ■ {a 2 h + {vi - V 2 )) \(t>{h)f fo{x - viti,vi) /o(t - V 2 ti,V 2 ), 


(4.19) 
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and, on the other hand (with /i 2 = h), 


^ ^ 0^02 J J '^'^1 J '^'^2 

^ ^ cri,(T 2 =±l 

f 1 - (Ji 1 + Ti ,(J2-(Ji 

0\V- Vi --- V2 —z- h 


2 2 2 
A+(cr2/i ■ {-h -Vi+ V 2 )) 4>{h) (j){-h - Vi + V 2 ) 
fo{x - viti, ni) /o(t - ^ 2 ^ 1 , ^ 2 ). 

Here we dehne the distribntion 


A+(a):=^ dse*“" =7rh(a)+ip.v. . 


(4.20) 


(4.21) 


Now, on both preceding formnlae, we readily observe the following important fact. The 
parity of (/>, and the symmetries h —h, ai — ti, cr 2 — ^ — cr 2 in (4.19), and h —h, 
CTi ——CTi, a 2 —^ —cr 2 , vi V 2 in (4.20), show that A+ may be replaced by nd everywhere. 
This will eventnally give, as shown next, the desired conservation of energy in the limiting 
U-U eqnation. 

There remains to actnally perform the snm ^ cria '2 in (4.19) and (4.20), in order to 

fl 21 (2t ll fl 2l 

identify the very valne of and T 2 ’ . For T 2 we make the following choice: 


Tl 

0-2 

h 

Vi 

V2 

1 

1 

v' — V 

V 

n* 

-1 

-1 

V — v' 

V 

n* 

1 

-1 

v' — V 

v' 

K 

-1 

1 

V — v' 

v' 

K 


This resnlts in the hnal expression: 

T 2 ^’‘^\x,v) J J dv' J dv'^ 5{v^ + V — v'^ — v') 

S + v^ - - n'^)^ \^{v' - v)\‘^ (/7* - ff*): 

where, with abnse of notation, we set the ’’transported qnantities” 
f = fo{x-vti,v), /* =/o(t f' = fo{x-v'ti,v'), /' =/o(a: 


(4.22) 




Notice that, by changing v' n', and nsing the conservation of momentnm, we may 
replace — n)p by — f*)!^ in (4.22). 

Besides, for make the following changes of variables: 


Tl 

T2 

h 

Vi 

V2 

1 

1 

V — v' 

n* 

V 

-1 

-1 

v' — V 

V 

n* 

1 

-1 

v' — V 

< 

v' 

-1 

1 

V — v' 

v' 

K 
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This results in the hnal expression: 

I n n 


^2 


{x,v) 


47^2 


(in* / dv' / (in* h(n* + n — n* — n') 


1 


5 ( -(n* +n^ 


n^ — n 


) ) ^(n' - n) ^(n' - n*) (/'/* - //*) 


(4.23) 


As a conclusion for the T 2 terms, we have eventually established the (desired) equality 
^a, 2 )^^^( 2 ,i) ^ j dv,w{v,v,\v\v:)iffi-ff,). (4.24) 

This ends up the analysis of the two-particle terms. 

5. Three-PARTICLE terms: permutations with a fixed element 


In this section we analyze W 3 g and V 3 g for the permutations tt with a hxed element. 
To simplify the notation we set 


W 3 g and V 3 g for tv = (1, 2, 3), 


and 


Wl, andV^,g, i = l,2,3, 

for the three permutations leaving i hxed. To state the result briehy, let us readily say 
that the factors 

^ 3 ,s: together with W| g, Vg^ g, 

give a vanishing contribution. Also, the sum 

>^3% + V|,g 

is shown to vanish asymptotically, while each of these two terms is 0{e~^) separately. 
Here, anti-symmetry will play a central role. Finally, the two terms 

wl,., vjg, 

do contribute to the limiting U-U equation through the cubic term. They build up the 
transition kernel (/)(p' — n)^ -|- — n*)^. The missing cross term 2(j){v' — v) (p{v' — n*) in 

IF(p, p*|p', p*) will come up in the next section. 

Let us show hrst that W 3 g and V 3 g are vanishing. From (3.18), scaling hi and /i 2 and 
summing on cri, (T 2 , we have 


^ 3 % = +4 


(^^ I I I *^^2 ^(£hi)^(£/l2) 


sm I sin ' {v + (ti - t 2 ){C - hi) 

MV. - (hi + h2),rj + ti^ - (tihi + t2h2)) Mhi.tihi) /o(/i2, ^2(^2) 


( 6 . 1 ) 
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due to the decay properties of /q. The same argument easily leads to = 0{e). 

We now pass to the computation of VV 3 g. This term is associated with the permutation 
TT = (1,3,2). Upon Fourier transforming in x the relation (4.2) for /J (with j = 3), and 
using the change of variables yi = {x 2 + xs)/2, 1/2 = {x 2 — xs)/s in the corresponding 
formula, we recover 


= £^/o(6,i7i) / dyi I dy 2 


+ f 2 / 2 ) p-j?3-(yi-f 2 / 2 ) 


172-173 , h 2 + 173 

yi - £—^—, -7/2 + —^— 


7 f . m-m ,172 + 173 
Jo It/I +£ —^—,T/2 + —-— 


Then, inserting (5.2) in the formula (3.18) relating the value of WJg, we arrive at 

aicr2 J dt2 j dhi j dh2 j dyi j dy2 (f){hi) (I){h2) 

7 hi + h2 , . f tihi+t2h2\ 

W-£-.^ + ^ 1 ^-^- 1 

~ / tl/ll — ^2^2 tihi -\-t2h2 \ 

fo[y, -J-,-- y.) 

7 f , tihi-t2h2 tihi+t2h2 , \ 

/o !/l +-J-.--+K . 


Changing variables 


we obtain 


hi + /l2 , , ti — t2 

k = -, h = 112 , s — -, 

£ £ 


YI Ticr 2 J ds J dh J dk j dyi J dy 2 (j){—h + sk) (l){h) 


.^^n■{-h+ek) ^-i^h-{r^+sh+ss{^-ky) ^-ik-y^-iy2{-h+e^) ^ 


tih kti + sh kti — sh 


\ 7 f tih kti + sh kti — sh \ 

T/2 1 /o ( T/i - ^ + ^- 1 -’-^-+h2 1 • 


We are now in position to identify the rigorous limit of W 3 g, using the assumed decay of 

(f) and /q. The argument are those used in the previous section: we do not repeat them 
here and in the sequel. Passing to the limit we get, eventually. 


wk(.(,y) = 


cricr 2 / ds dh dk dyi / dy 2 \4>{h)y 


<ti,(T2=±1 




~ / tih kti — sh 

h[yi + —.^— 


- T/2 1 /o ( T/1- 


tih kti — sh 
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Last, we go back to the {x,v) variables, computing the inverse Fourier transform of the 
above term. This gives 


Wl{x,v) = -TT ^ aia2 / dvi / dv2 / dvs \(f){v2 -1^3)!' 
cri,cr 2 =±l J J J 


d[[v2- 1 ^ 3 ] ■ 


-V + V2 


1 + CTi 


+ Vs 


I- ai 


2 '"2 
/o(t - Uiti, Ui) /o(t - V 2 tl, V 2 ) fo{x - Vsti^Vs). 


d (v - Vi + {vs - V 2 ) 


(5.6) 


This is our hnal expression of VV 3 . It will be interpreted later in terms of the v, n*, v', n' 
variables of the U-U equation. 

In a similar fashion we compute V|g and its limit V|. We write 


2/1 - £ 


vi - m 


, - 2/2 + 




hi +h3 


V 


4 2 


(5.7) 


We insert this expression into (3.19), and perform the change of variables h = h 2 , k = 
{hi — /i 2 )/£, and s = {ti — t 2 )/s. Passing to the limit £ —0 at once, gives the asymptotic 
value 

V|(^,h) = - 77 ^ ^ 1^2 I ds I dh I dk I dyi j dy 2 \${h)\^ 


(27r) 


<ti,(T2=±1 


2/1 + 


tih rj + ^ti — kti — sh 

2 


- 2/2 ) /o I 2/1 - 


tih T] + — kti — sh 


+ 2/2 


fo{k,tik + sh) e 'v g--- t-'—/ g 

whose inverse Fourier transform is 

V|(T,n) = -7r ^ aicr2 / dui / (in2 / ^ 1^3 |^(i^i-w)h 

<ti,<t2=±i j j j 


(5.8) 


5 ( [ui - Us] 


f 2 - 'I’l 


1 + 0-2 1 - 0-2 
- 1^3- 


2 2 
fo{x - Viti,Vi)fo{x - V2tl,V2) foix - Vsti^Vs). 


, 1-0-1 1 + Tl 

o\v- Vi —-- Vs --- 


(5.9) 


Before coming to the computation of the other W^’s and V^’s, let ns now identify the 
link between the obtained values of VV 3 , Vf, and the U-U equation. The following changes 
of variables in VV 3 


Tl 

T2 

1^1 

V2 

Vs 

1 

1 

V 

K 

n* 

-1 

-1 

V 

n* 

< 

1 

-1 

v’ 

< 


-1 

1 

v' 

n* 

K 
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yields the explicit value 


yV^{x,v) = 2 tt j dv^ / dv / dv[ \(p{v —v 


For Vq we set 


Tl 

0-2 

1 

1 

-1 

-1 

1 

-1 

-1 

1 


h(u + u* -v' -v'^) 5 [ -iy^ + vi -v'^ - < ) [/*/7* -//*/*] 


( 6 . 10 ) 


and obtain the hnal 


Vl{x,v) =271 / dv^ / dv' / 7 (n-n')p 


s(v+v, -v' - v',)s [//'/; - /;//.]■ 


( 8 . 11 ) 


Last, using the symmetry v' n' (note that v' — v — —{y'^ — v^) and v' — v^ = 'i’ — 'I’D we 
hnally conclude after some computations that: 


{Wl +Vl){x,v) ^ Ti J dv^ J dv' j dv'^ [(/ +/*)/ 7 *-(/' +/*)//*)] 

S{v + v^-v'- ifl)6 + ''j^ - ilHV - ^')l^ + \y' - ^'*)l^). 


(5.12) 


This is the desired cubic term in the U-U equation, up to the fact that we only recover 
here part of the total cross-section W{v,v^\v',v'^) — [(f){v — v') -|- 4>(y — The missing 

cross term will come up in the next section. 

We now show that all other terms associated to permutations with a hxed element, 
namely W| g, V 3 g, and Wf ^ -|- Vf ^, give a vanishing contribution in the limit £ —0. 

We begin with Wfg. Inserting (5.7) into (3.18) and changing variables k = 

/i 2 = h, we readily obtain: 

aia '2 [ dt 2 [ dk [ dh [ dyi [ dy 2 4>{h) (f){ek) 


^ Jo I 

Q-^{o-ivke+a2h-{r]+{ti-t2)y-k)) ^iy2-h g-*(yi+e^)-(?-^) 

~ / £ r 2 1 1 

/oh/1-7 V + Ch-tik - t2h ;-{v + Ch - tik) - y2 

\ 4 L £ J 2 

~ f £ r 2 1 1 

foiyi + ^ Tj + - tik --t 2 h ■,-{rj + ^ti - tik) + y 2 

fo{k,tik). 


(5.13) 
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Summing on ai and <72 allows to compute the limit which is: 
4 




(2t) 


dh dk dh dpi / dy^ ,#(h) .^(0) ^ 


sin Qh ■ (rj + {ti - t 2 )(^ - k)^ fo (^Vi + ^^ 2 / 1 ; ^(ry + ^ti - tik) - y 2 

fo (yi - ^t2h; ^{v + ^ti -tik) + 1 / 2 ^ fo{k,tik). 


(5.14) 


Note that the product of the two /o’s is invariant under the change of variables h —h, 
1/2 —y 2 : as well as the oscillatory factor Therefore, using the parity of (/>, it follows 

that Hence ^ 

Wi = 0 . 

The term V 3 is treated in the analogous way. We insert (5.2) into (3.19) and make the 
change of variables h = /i 2 , k — e~^hi, obtaining after summation over ui and cr 2 : 


^^2 J J J j ^y^ ^{h)^{ek) 


^-iyi-k iy 2 -{h-e^) 


f £ 

2 ■ 

(7/1-7 

tik - t2h 

V 4 

e 


sin (sin ^ ■ k{ti - t 2 ) 


1 

-! 

2 


f £ 

2 ■ 

(7/1 + 7 

tik - t2h 

V 4 

e 


; -^tik + y2 


foiC - k,r] + ti{^ - k)). 
This term clearly goes to 
4 


Va* «,<?) = 


(27r; 


6 


dt2 dk dh dyi / dy2 (j){h)(j){0) 


>0 


^-ryi-k ^dy 2 -h sin (^h ■ k{ti - t 2 ) 


1 ,1, 

yi + -t 2 h-, -tik-y 2 


fo (^yi 


^t2h-,^tik + y2] fo{^- k^y + tii^- k)). 


(5.15) 


(5.16) 


Hence V 3 = 0 for the same reason as before. 

To end up this paragraph, let us last prove that the sum Wf g + Vf ^ vanishes asymp¬ 
totically. First, we write 


.^,e(6,6, 6, 171 , 172 , 173 ) = s'*/o( 6 , 173 ) / C^77i / dy 2 e 


-i^i-(yi+§y2) „-i?2-(yi-f2/2) 


fo[yi-e 


m - m 


,-7/2 + 


m +172 


/o 7/1 + £ 


m - m 


,7/2 + 


m +172 


(5.17) 
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We insert this formmula in (3.18), and perform the change of variables hi = h, h 2 = sk. 
In this way we recover, npon compnting the snm X] ^ 1 ^ 2 , 

^^2 j dk j dh j dyi j dy 2 ^{h)^{ek) 


QW2-h Q-i{yi + ^y2)-{i-k) 


sin -?7 ■ h sin — ■ (?7 + {ti - ' Kh - ^ 2 ) 


C ^ J. 

2 / 1-7 + - tih ;-{h + Cti-t 2 k)-y 2 


2/1 + 7 ?7 + ^ti-H 2 - hh ;-{v + Ch - t 2 k) + y 2 


fo{k,t2k). 

(5.18) 

Similarly, nsing (5.17), (3.19), and performing again the change of variables hi = h, 
/i 2 — sk, we obtain: 


ViA(,rt]=- 


--1 


dt2 dk dh dyi / dy2 (f>{h) (f){ek) 


Qiy2-h g-*§2/2• (?+fc) 


sin ( -?7 ■ h ) sin ( -k ■ h{ti - 12 ) 


/ £ 

2 ■ 

2/1 + 7 

V + ^ti + kt2 - tih 

V 4 

e 

(yi - J 

y + Cti + kt2 - -tih 

V 4 

e 


(5.19) 


fo{k,t2k). 

Hence, both terms Wfg and Vfg are 0{e~^). However we have the following expansions: 

^^3,£ — ^0 + £^i + 0{e‘^), eVfg = Bq + sBi + 0{e‘^), 

and it is easy to realize that Aq = —Bq. Moreover, after some straightforward calcnlation, 
we obtain at the next order: 


- —1 /"ii 


Ai+Bi = 


dt 2 dk dh dyi / (ii /2 </>(h) (/>( 0 )e*^^'^ e sin i-r] ■ h 


1 , 1 


fo{k, t 2 k)fo yi - -tih; -{y + ^ti - t 2 k) + 1/2 /o 2 /i + 7 ^ 1 ^; 7 ;iv + CA - t 2 k) - 1/2 


2 ’ 2 


■ {rj + (ti - 22)0 cos ( -k ■ h{ti - t 2 ) + 


1 , 1 


2 ’ 2 


sin -/c ■ h(ti - 22 ) 7 ;t 2 k-Vxiog 


foUvi - \tih-, \{ri + ^ti - t2k) + 1/2) 

7 o(( 2 /i + \Ah] i(?7 + ^ti - t2k) - 1/2) 


-iy2-k 


(5.20) 
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Let US now exchange h —h and y 2 — 2 / 2 - The term in braces is invariant because the 
log term changes its sign. All the other terms are invariant but sin (?7 ■ h/ 2 ), which changes 
sign. Therefore Ai + Bi = —{Ai + Si), hence Ai + Bi = 0. This shows that Wf ^ + Vf g 
vanishes in the limit £ — 0 . 


6. Three-PARTICLE terms: cyclic permutations 

We still have to evaluate WJg, ^ ^3 s '^3 with tt = (2,3,1) and 7 r“^ = 

(3,1,2). 

We hrst observe, for later convenience, the two relations 
fs (6,6,6;hi,h2,h3) = [ dxi [ dx 2 [ dT 3 


, , Xi + X2 e , . Xi - X2 Vi + r]2 

fo I -^-+ 4 W “ ^ 2 );-^^- 

'X 2 + X 3 e . X 2 - X 3 ri 2 + rj 3 

fo 1 -^-+ |(h 2 - 173 );-^^- 

'X3 + X1 e . X3 - xi r]3 + rji 

fo 1 -^ 4 (h 3 - hi); —^^— 


/:r(6,6,6;hi,h2,h3) = / dxi / dx2 / dx3e 


^ Sfe=i 


( 6 . 1 ) 


,.X2+Xie . X2 - Xi rj2 + m 

fo I —^— + ^(h 2 - hi); —^^— 

, 'X3+X2 e . X3-X2 Vs+ rj2 

fo 1 -^- + |(h 3 - h 2 ); —^^— 


xi+xs e . X1-X3 rji+rjs 

^—+ j(-?i - ^—+ 


( 6 , 2 ) 


Armed with these expressions, we begin with the computation of VV 3 g and WJg. To do 
so, we insert (6.1) and (6.2) in the general formula (3.18) relating the value of the WJg’s. 
In the so-obtained formulae, we also change variables, hi —hi, /i 2 —^ —/i 2 for 7 r“^, and 
ai —— Ti, (72 —^ —(72 for TT. With this new set of variables, the ^’s and ry’s involved in 
(3.18) are 

^ ^ . hi h2 . hi /i2 

-, 6 = t—, 6 = t—, 

£ £ £ 


I tihi-\-t 2 h 2 tihi 72^2 

hl=h + ?fl± -, h2 = T-, h3 = T-: 


(6.3) 


for TT ^ and tt respectively. Also, the phases appearing in (3.18) are: 


= y/ii ■ h + y ^2 ( h + (^1 - ^ 2 ) 




hi 


Sn = y hi ■ Tj + y h 2 ( h + (^1 - ^ 2 ) 




(6.4) 
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All this gives in (3.18), the two values 

^3,£ =- cria 2 j dt2 j dhi j dh2 j dxi j dx2 j dx3 

~ fxi+X2 2 tihi+t 2 h 2 £„ Xi - X2 ^2^2 f] 

* +-4-+ 4’'’^— +^+2 


/o 


X2 -h X3 tihi — t2h2 X2 — X3 tihi -\-t2h2 

2 4 ’ £ 2£ 


, I'xs + Xi tihi+2t2h2 s_ X3-X1 , tihi , r] 

-4-+ 


(f){hi) (f){h 2 ) e 


-4X1 ■(g+ A . 


(6.5) 


T 1 CT 2 J dt 2 J dhi J dh 2 j dxi J dx 2 J dx 3 

~ fx 2 + xi 2 tihi+t 2 h 2 e_ X 2 - xi ^ 2/12 V 

+- 4 - 4^’^—-^+2 

J [X 3 +X 2 tihi-t 2 h 2 X 3 -X 2 tihi+t 2 h 2 

- 4 -’ +- 2 i- 

~ f xi + X 3 tihi+ 2 t 2 h 2 e_ xi - X 3 tihi f) 

* ^- 4 — + 4’'-^—ir + 2 


(j){hi) (j){h 2 ) e 


-4xi-(g- 


( 6 . 6 ) 


where we use the notation rj = rj + ^ti. 

Now, we perform the following natural change of variables: 


for TT ^ 
for TT : 


1 , 1 , 

T2 = Ti + -t2h2 - eyi, X3^Xi- -ti/li + £1/3, 

1 , 1 , 

X2 = XI + -t2h2 + eyi, X3 = xi - -tihi - £1/3. 


(6.7) 


In both cases xi is unchanged. This hnally gives the two relations 

- — 1 rti 


7^3,5'(^, 17 ) ^ 1^2 j dt2 J dhi J dh2 J dxi J dyi J dy3 


criO-2=±l 


~ f tihi + t2h2 , £ 
Jo I Tl H--h 


/o Ti - 


a 

2 ' 2 12 
tihi — t2h2 £ 


yi 


;yi + 


V 


- 2 ( 2/1 “ 2 / 3 ); -yi - ys 


7 f tihi + t 2 h 2 £ 

/o (^4,-^-- 


'd . 

2 


V 

;i/3+ 2 


(f){hi) 4>{h2) e 


^-iyi-hi ^iy3-h2 ^-^hi-h2(ti-t2) 


( 6 . 8 ) 
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dys 


fo{xi + 


cri(T2=±l 

tlhi + ^2^2 £ 

2 2 


y 

2-y\ 


, ^ 
;2/i+2 


~ / ti/ii -^2^12 X 

Jo I Ti-^-+ -{yi - 2/3); -yi - ys 


(6.9) 


2 L2 


- 2/3 


;2/3 + 


V 


^{hi) ^(/i2) 


Next, we come to the computation of V 3 ^ and V 3 ^ . We insert (6.1) and (6.2) in (3.19). 
We also change hi —hi, ct 2 —^ —cr 2 for and ^2 —^ —^ 2 , cr 2 —^ —cr 2 for tt. This gives 
the two identities 


vj: (e,i7) = - 


--7 


( 271 ) 


E 

(Ticr2=±l 


fti 


cria2 


(it 5 


dhi / dh2 / (i^i / (ia;^ 


(ix? 


/o 

7o 


/a:i+a:2 2tihi+t2h2 e_ xi - X2 ^2/12 17 

+-4-* 4 ’'’- TT + 2 

fx3 + xi tihi-t2h2 e _ X3 - Xi tihi+t2h2 fi\ 

[— -4-*-27-* 2 ) 

^a:2+T3 tihi + 2t2h2 X2 - X3 tihi 
^ ^ ^ 


V 2 

7(hi) (j){h 2 ) e 


4 


( 6 . 10 ) 


V 3 " (e,i7)=- 


--7 


( 277 ) 


E 

cri(T2=±l 


fil 


Ti(72 


dt2 / dhi / dh2 / (ixi / (ia;2 / (iT3 


~ /a :2 + Ti 2 tihi+t 2 h 2 e_ X 2 - xi t 2 h 2 y 

+- 4 -+ ^+2 

T 1 +T 3 tihi-t 2 h 2 e_ xi - X 3 tihi+t 2 h 2 f) 

— - 4 -+ 4 ’'-^- 27 - *2 

~/a;3 + a :2 tihi+ 2 ^ 2/12 X 3 -X 2 tihi 

^^[~~2 - 4 -+ ^ 


^(hi)^(h2)e 

Here the phases are: 


g-ja:2-E±l:2 


£ e £ e 


‘5^-1 = ^hi ■ y + ^h2 ■ —(ti - t2), 
2 2 e 

^ CTi (72 hi 

= -;r"i ■ h - -Trh2 ■ —(ti - 12). 

2 2 e 


( 6 . 11 ) 


( 6 . 12 ) 
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We make the following natural change of variable 

1 1 / \ 
for TT~^ : xi = x'l + -^ 2 / 12 , X 2 == - eyi, = x\ - -jtihi - £ (^i/i + 1/3 + ^ j , 

for TT : Xi^x[ + -t 2 h 2 , X 2 ^x[+eyi, X 3 ^ x[ --tihi + e {yi + ya +-J . 

(6.13) 

With this change of variables, we eventually obtain 

o'icr 2 j dt 2 J dhi J dh 2 J dx'i J dyi J dys 

cri(T2 = il 


fo(xi + 


, tihi + t 2 h 2 e r ?7 


e r ?7 

2 [2 ’ 


, V 

!/i + 5 


/ tihi — t2h2 £/- , , X 

-^- ^{v + yi + 2/3); -yi - 2/3 


~ / , tihi +^2^2 £ r ?7 

/oO,-^- 


£ r?? ^ 1 

- 2 [ 2 ^ ^ 2/3J; ys + 


${hi) 5(/i2) e^ys-h2 


(6.14) 


as well as 


Kei^,V) = 


(J 1 CT 2 J dt 2 J dhi J dh 2 J dx[ J dyi J dys 

(Ticr2=±l ® 

hi +t 2 h 2 e \fj 1 fj\ 

2 2 [2 2 j 


r / / , +t 2 h 2 e \r] ] r]\ 

k[xi + -^-2 [2“H’^^+2j 

7o (x[ - ^ ^ 2/12 ^ ^ 


fo{xi- 


, tihi + t 2 h 2 ,£ Ih 


'V rx 1 

- + 2 yi + ysj ; 2/3 + 


(/)(/ii)(/)(/i2)e 


iyz-h2 „-^hi-h2{ti-t2) „i^-{v+^hl-t2)) O-Stt 


(6.15) 


Let us come to the computation of the limit of the above four terms WL, , WL^, VL, , 


VJg. The phases carried by these terms are respectively: 


aihi+a2h2 , , .r, ,x . , , , 1 - ^2 ti -12 , , 

- - y + —h2 ■ 4 (fi - ^2) - Ti ■ 4 - ■ hi + 1/3 ■ /i2--- hi ■ h2, 

2 2 2 £ 

aihi+a2h2 , cr2, , ,x ^ . 1 - ^2 ti - ^2 , , 

- - y + —h2 ■ 4 (fi - ^2) - Ti ■ 4 - yi ■ hi + 1/3 ■ h2 H---hi ■ h2, 

2 2 2 £ 

aihi + h2 1 , . , , l + T2ti-t2, , 

-- V + 0^2 ■ 4 (fi - ^2) - ■ 4 - I/I ■ hi + 1/3 ■ h2 H-r- hi ■ h2, 

2 2 2 £ 

aihi + h2 , Ij, ^ I, , i, l + T2ti-t2j, J, 

-- V + oh2 ■ 4 (fi - ^2) - ■ 4 - yi ■ hi + 1/3 ■ h2-r-hi ■ h2. 

2 2 2 £ 
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Denoting by Wg ^ yV^’^, Vg ^ the eight terms relative to the valnes of a2 = ±1, 

we realize that Wg ^ Wg Vg ^ Vg g~, have only slowly varying phases, so that they 

are individnally 0{e~^). However, setting 

= ^0 + + 0{e^), eWg’g^ = Bo + Bie + 0{e'^), 

__ _ (6.16) 

sVg^g ’ = Cq + Ci£ + 0(£^), = Dq + D\£ + 0(£^), 

an easy hrst order Taylor expansion gives Hq + C'o = Bq + Dq = 0 and Hi +i?i = Ci+Di = 
0. Hence 

lim (>vj;‘’++wjy+vj;'--+ %■-) = o, (e.i?) 

For the other terms, which carry a rapidly oscillating phases, it is natnral to rescale time, 
setting s — £~^{ti — t 2 )- Then, an easy compntation shows 

1™, (tvjy + =■ tvj. (6.18) 


where 


f dh 2 ^{hi) ^{h 2 ) 


^—ihi-h2S _j_ ^ih\-h2S 


dxi / dyi / dyg 


/ 


/o( 


hi + /i2 , ?7 

+fi—^—;2/i + - 


r / , hi + ^2 

Jo 1 Ti -ti -^-;y3 + 2 


2 y ■^° ( 
\ 

V 


Xi — ti 


hi - h2 


■; -2/1 - 2/3 


(6.19) 


Finally, taking the inverse Fonrier transform of this term, we obtain: 


Wg (T,n) =27r j dvi J dv 2 J dvs (t>{v 2 - vi) (f){vs - V 2 ) 

[h(n + ^2 - ni - V3) - d{v - na)] d{{v2 - ni) ■ (ns - ^2)) 
fo{x - nHi, vi)fo{x - V2ti,V2)foix - ngti, ng) 

=27r J dn* J dn' J dv' ^{v' - n*) ^{v' - v) {/*/'/* - //*/'} 

S{v + n* - n' - n') 6 


v'^-vT, 


( 6 . 20 ) 


In the similar way we compnte the limit 


lim 

e^O 


V. 


3,£ 


+ v. 


71,+ 

3,e 


= V. 


3 5 
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whose inverse Fourier transform admits the value 


V'^{x,v) =27r j dvi J dv2 J dvs (I){v2 - vi) (l){vs - V 2 ) 

[h(n - V2) - S{v - ui)] d((v2 - vi) ■ (V3 - V2)) 

fo{x - Viti, Vi)fQ{x - V2ti,V2)fo{x - Vsti, V3) 

=27r J dv^ J dvi J dv'^{v'- - v){f f fi - f f^fi} 

d{v + v^ -v' -v'JS 


( 6 . 21 ) 


There remains to sum up the contributions of the terms WJ and VJ. It gives, after 
some computations using the exchange of variables v' ^ vi, the missing cross term 

(W 3 + VJ)(a:, v) = 27r J dv^ J dv' J dvi [if + - if' + /*)//*)] 

d{v + v^-v' - vi) 5 +vl- v'‘^ - ^{v' - v) 4){v’ - n*). 

This completes the proof of the theorem. 


7. Concluding Remarks 


It is well known that other possible scaling lead to kinetic equations as well. The most 
important is the low-density limit (or Boltzmann-Grad limit): it is the regime in which 
classical rarehed gases are described by the usual Boltzmann equation. 

In our grand-canonical formalism it can be introduced in the following way. We do not 
rescale cf) which is 0(1). On the other hand the rarefaction hypothesis is given by the 
condition e^{N) = 0(1). This means that (see (2.24)) 

f^{x,v) = 0{e). (7.1) 

Under this scaling, the hierarchy becomes (see (2.18)) 


dj’ + Y.V,- = R'/' + 

k=l 

Rescaling the correlation functions by dehning: 

we arrive at the hierarchy 


e 

3 + 1 ' 


dj‘ + ■ V,J| = iqyj + 

k=l 


(7.2) 


(7.3) 


(7.4) 
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with a fixed initial datum of 0(1). It is now clear that the terms CC are vanishing in 
the limit £ ^ 0 and the statistical correlations are lost. On the other hand many terms 
of the type CT.. .T are finite in the limit. It turns out that the sum of these terms 
lead to the Born series expansion of the cross section. The underlying series actually 
converges provided the potential (f) is small. This task is performed in the case of the 
Maxell-Boltzmann statistic in Ref. [3] by the authors. Here, a difficult point lies in the 
identification of the cross section as the Born series expansion of quantum scattering, a 
task which is achieved using an original identity derived in [6]. 

Another comment is in order. The U-U equation has been partially derived whenever 
/o is the Wigner transform of a one-particle quasi-free state. As shown in Appendix, a 
sufficient condition for the explicit construction of such a state is a small value of the 
activity 2 . On the other hand the U-U equation for Bosons makes sense also for more 
general initial conditions describing states with large activity. It seems very interesting to 
understand whether the U-U dynamics of such states make sense from a physical point of 
view and whether it can describe dynamical condensation phenomena. 

Appendix: Quasi-free states for Bosons 


Let r be a one-particle state i.e. a self-adjoint positive operator whose kernel is denoted 
by r{x,y). We want to construct a state which is compatible with the B-E statistics and 
with a given average particle number. 

Let an be a n-particle completely symmetric state given by 


a. 


^ ^ {parity-K^n))- 


TreVn 


The state 


where 


a" = 


1 °° 

^ ' n=0 




n>0 


is a normalized state for Bosons and 


d 


(TV) = Tr [a^iV] = logS. 

dz 

We now compute the partition function S(z) (see [9]): 

= [ dxi...dXn r{xi,x^(i)) . . .r{Xn,X^(n))- 


1>0 


(A.l) 


(A.2) 


(A3) 


(A.4) 


(A.5) 




Given tt, let ai,..., an, be non negative integers, aj denoting the number of cycles of 
length j in TT. Clearly 



ON THE WEAK-COUPLING LIMIT FOR BOSONS AND FERMIONS 


29 


Given the sequence ai,..., 


dxi ...dxn r{xi,x^(^i ))... r(a;„, = JJ (^Trr 

j=i 


The number of permutations associated to a given sequence cti,..., is 

n 

nlTT—r—. 

J. J. rv d 1 “i 


a-jlT 
i=i ^ 


Hence 


E n 

n>0 = l 

'^jaj=n 


Trrdj 


Tr 


EEs EH,, 

'^jk=n 


= exp 


E 

j>i 


Trr^ z^ 


J 


(A7) 


(H.8) 


(H.9) 


In the second equality s denotes the number of actual cycles in each permutation and 
ji,... ,js are the lenghts of the cycles. The last sum is convergent for ^ sufficiently small 
(away from Bose condensation region). Then, by (A.4) and (A.9) 


(N) = ^ Tr {rzy — z + o(z), 
j>i 


(A.IO) 


for z small. 

The RDM’s according to (2.5) are: 

1 („ -u D’ f 

= SS) ^J (. 4 , 11 ) 

Therefore we have: 

Pji^j^ y?) ~ N ^ ^ 7 _ PI ^ ^ + ' ' / dZn r(xij . . .r(Zm CTr{n))i (A.12) 

where ^ = {Yj^Zn-j). Hence 

Pj{Xj,Yj) = ^ yy yy f ^ j ^ dZsr{zi,z,,/(^iy...r{zs,ZT,/(^s)) 

^ n>j ^ s=0 ^ 2 Tr'ePs 

(n-i-s)!^ ^ Wr^yxi,y^^i)) 

TvE'Pj ki,...,Kj'>l i=l 

J2kj = n-s 


(. 4 . 13 ) 
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with r^{x,y) the kernel of 
Since 

n- j\ {n - j - s)l _ 1 
s 


{n-j)l s!’ 


^sn — s T,ki 

— z — z z — z z 


we obtain 


3 oo 


' s>0 *■ TrePi £=1 ke = l 


Dehning the one-particle operator 


Tz = z'~r'~ — 

k>i 


^kk 


zr 


1 — zr ’ 


for 2 snfhciently small, we arrive at 

'y ^ z{,^i-)y'K{i)) • ■ • z{p£j•)y-K{j))•) 

■neVj 


(A.14) 


(A.15) 


(A16) 


(A.17) 


that is the characterization of the qnasi-free state in terms of RDM. 
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